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Abstract 



> 
in 

CO 

CN ' We apply the Ginzburg-Landau theory to the colour superconduct- 

ing phase of a lump of dense quark matter. We calculate the surface 
, energy of a domain wall separating the normal phase from the super 

phase with the bulk equilibrium maintained by a critical external mag- 
^ netic field. Because of the symmetry of the problem, we are able to 

simplify the Ginzburg-Landau equations and express them in terms of 
two components of the di-quark condensate and one component of the 
gauge potential. The equations also contain two dimensionless parame- 
ters: the Ginzburg-Landau parameter k and p. The main result of this 
paper is a set of inequalities obeyed by the critical value of the Ginzburg- 
Landau parameter-the value of k for which the surface energy changes 
sign-and its derivative with respect to p. In addition we prove a num- 
ber of inequalities of the functional dependence of the surface energy on 
the parameters of the problem and obtain a numerical solution of the 
Ginzburg-Landau equations. Finally a criterion for the types of colour 
superconductivity (type I or type II) is established in the weak coupling 
approximation. 
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1. Introduction. 



The Ginzburg- Landau theory and the surface energy 



The Ginzburg-Landau theory provides a powerful tool for exploring systems with in- 
homogeneous order parameters near the critical temperature. For instance, its application 
to the surface energy of a normal-superconducting interface led to the discovery of the two 
types of superconductors. When an external magnetic field whose magnitude is less than a 
critical value H c , is applied to type I superconductors, the field is expelled from the interior 
of the sample (the Meissner effect) . As the magnitude of the magnetic field increases above 
H c , superconductivity disappears completely and the normal phase is restored. In the case 
of type II superconductors, the sample behaves in a similar manner as long as the external 
magnetic field remains below a lower critical magnitude, H c \. A qualitatively different 
behaviour emerges when the applied magnetic field lies in the range H c \ < H c < H C 2, 
where it becomes energetically favourable for the magnetic field to penetrate the sample 
in the form of quantized flux lines called vortices. Beyond the upper critical magnitude of 
the field, the vortices are too dense to maintain the condensate, and the normal phase is 
once again restored. 



The criterion that determines whether a superconductor is of type I or type II is the 
Ginzburg-Landau parameter. It is defined as the ratio of the penetration depth of the 
magnetic field 5 over the coherence length £, the distance over which changes in the order 
parameter occur, i. e. k = |. These characteristic lengths imply that fluctuations in the 

magnitude of the condensate decay as e _v ^f while the magnetic field inside the super- 
conductor falls off as e~% . In their original paper, Ginzburg and Landau [1] analyzed the 
energy of the interface between a normal and a superconducting phase, kept in equilibrium 
in the bulk by an external magnetic field at the critical value. They found analytically 
that the surface energy vanishes at 

k = k c = —j= — 0.707. (1.1) 
v2 

The physical meaning of this critical value was clarified further by Abrikosov [2] . It repre- 
sents the demarcation line between type I (k < -4=) and type II (k > -4=) superconductors. 



In the present paper we shall carry out a similar analysis of the Ginzburg-Landau 
theory for colour superconductors. Colour superconductivity is essentially the quark analog 
of BCS superconductivity [3], [4], [5], [6]. Because of attractive quark-quark interactions in 
QCD, the Fermi sphere of quarks becomes unstable against the formation of Cooper pairs 
and the system becomes superconducting at sufficiently low temperatures. In this paper we 
shall consider a domain wall separating the normal and colour superconducting phases in 
equilibrium under the influence of an external magnetic field and calculate the surface free 
energy per unit area. The Ginzburg-Landau equations for this system are considerably 
more complicated than the analogous equations for an ordinary superconductor [7], [8], 
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[9], [10], [11], because of the non-abelian nature of the problem. By using symmetry 
arguments, we were able to simplify the problem of the domain wall and express the G-L 
equations in terms of one gauge potential and two components of the di-quark condensate. 
When expressed in terms of dimensionless quantities, the G-L equations contains two 
dimensionless parameters: the usual Ginzburg-Landau parameter k and another parameter 
p. Weak coupling calculations fix the value of p to be equal to — |. Although we were 
unable to derive an analytical result analogous to (1.1) we proved the following inequalities 
for the critical G-L parameter and its derivative with respect to p: 

* w £ 7T t £ °" (L2) 

We also derived a number of inequalities of the functional dependence of the surface energy 
on the parameters of the problem and obtained a numerical solution for K c (p). For weak 
coupling we found that 

k c ~ 0.589, (1.3) 

in contrast to (1.1). 



In section 2 of this paper, we shall review the Ginzburg-Landau theory of colour 
superconductivity. The analytical and numerical solution of the domain wall problem will 
be presented in section 3. In the final section of the paper we shall summarize our results 
and in the appendix we shall discuss the validity of the Ginzburg-Landau theory in the 
presence of the fluctuations of the gauge field. 



2. The Ginzburg-Landau Theory of Colour Superconductivity. 



The symmetry group of QCD in the chiral limit is 

SU(3) C x SU(3) fR x 5^(3)^ x U(1) B , (2.1), 

where the subscript c denotes colour, the subscript /r(/l) the right (left) -hand flavour and 
B the baryon number. The electromagnetic gauge group U(l) em is not a separate symme- 
try group, but a subgroup of SU(3)f R x SU(3)f L . The dominant pairing channel consists 
of two quarks of the same helicity and the corresponding order parameters in a colour 
superconducting quark matter will be denoted by and ^>l- These condensates trans- 
form into each other under space inversions. Neglecting the parity violating processess, 
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the Ginzburg-Landau free energy functional which is consistent with the symmetry group 
(2.1) reads 

d 3 f{-F^Fl + -(V x A) 2 + -Tr[(D9 R )i(D9 R ) + (D^ L )^(D^ L )} 

+ ^aTr[4** + + \bTr[(^ R ) 2 + (tfj. * L f] 

+ ±b'[(Tr*{* R ) 2 + c&*t^) 2 i 

+ ^cTr(^ R )Tr(^ L )}} 

(2.2) 

where the gauge covariant derivative of the di-quark condensate \& reads 

{D* R{L) )l% =V(v^ (i))? - - ^ c '(*^))^; 2 - igA^\* R{L) )X h 
-ie(q h + qh )A^ R{L) r f [% 

where A = A l T l denotes the classical vector potential of the SU(3) colour gauge field, A 
the electromagnetic field, and eqj is the electric charge of the f-th flavour quark. The 
SU (3) generator T l is in its fundamental representation. 

At weak coupling, the parameters are calculable from either the perturbative one- 
gluon exchange interaction of QCD or from the Nambu-Jona-Lasinio effective action. Both 
approaches lead to the same expression [7], [8]: 

-k B T c (T — T c ), 



7C(3) 

576tt 4 fk B T c \2 (2.4) 



b 7C(3)V p 
b' = c = 0, 

where fx is the chemical potential and T c the transition temperature. 

Since in this paper we address the domain wall problem- the interface between a normal 
and a superconducting phase in an external magnetic field- the boundary conditions select 
the even parity sector of (2.2) 

V R = y L = y. (2.5) 



The one-gluon exchange process that dominates the di-quark interaction at ultrahigh 
chemical potential is attractive for quarks within the colour antisymmetric channel. Thus 
the Cooper pairs realise the colour antisymmetric representation to the leading order of 
the QCD running coupling constant at T = [12], [13], or to the leading order of (1 — ^r) 
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near T c [8]. Assuming such a pairing pattern persists even at moderately high chemical 
potential, the di-quark condensate \1/ can be expressed in terms of a 3 x 3 complex matrix, 

= ^^hhf^f (2-6) 
The gauge covariant derivative of $ can be written as 

/ 

-> -t -±cc , 

(D$) c f = - igA $ c f -ieQ f A$ c f (2.7) 

where 

A = A l f l , <$> = <$> + <$> l f l (2.8) 

with T l = —T l * being the generator of the 3 representation, Qf = + qf 2 and //1/2 
represent a cyclic permutation of 1, 2,3. Arranging the flavour index in the conventional 
order of u, d, s we find the diagonal electric charge matrix 

Q = - dia 9i\-\-\) = -^f\ (2.9) 

Here we have adapted an expression of T l that differs from the standard one by a cyclic 
permutation of rows (columns), in which T 8 = -^K^diag(2 1 — 1, — 1). The equation (2.7) 
takes the matrix form 

£)<$> = V$-igA$ + i-?=eA$f 8 , (2.10) 

v3 

and the Ginzburg-Landau free energy (2.2) becomes 

/r 1 1 
d 3 f -F^F-j + -(V x A) 2 + 4tr(z3$)t(l3$) +4atr$t$ + 6 1 tr($t$) 2 + 5 2 ( tr $t$) 

(2.11) 

where the parameters b\ and 62 are related to the parameters of equation (2.2) in the 
following manner 

&! = &, b 2 = b + 8b' + 8c. (2.12) 



First, let us review the case with a homogeneous condensate, i.e. A = A = V$ = 0, 
for which the Ginzburg-Landau free energy (2.11) becomes 



= J d 3 f 4atr$t$ + b ltr ($t$)2 + b2 ( tr $t$); 



(2.13) 



We distinguish the following three regions of the parameter space b\ — 62, following the 
treatment of [7] 



1) bi > and b\ +3&2 > 0: The minimum free energy corresponds to the colour-flavour 
locked condensate [4], 

$ = (j)oU (2.14) 
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where 

*- = v /I ra (2 - 16) 

and U is a unitary matrix. ^ Consequently, we find that 

r min = -n-r^^r-. (2.16) 

2) bi < but 61+62 > 0: In this case the colour-flavour locked condensate (2.14) becomes 
a saddle point of the free energy, which is neverthless bounded from below. The minimum 
corresponds to an isoscalar condensate, given by 

$ = diag(^ e ia ,0,0). (2.17) 
where <j>' = yj — b ^ b2 • We find then that 

4a 2 

r min = -0 . (2.18) 

bi + 3o 2 

3) For bi and 62 outside the region specified by 1) and 2), the free energy is no longer 
bounded from below. Higher powers of the order parameter have to be included and the 
superconducting transition becomes first order. Since case 2) is mathematically identical 
to a metallic superconductor, we shall not address it in this paper. 



In order to identify the various characteristic lengths of the system, we consider fluc- 
tuations about the homogeneous condensate (2.14). We parametrize the order parameter 

.-^(^H^***) (2.19) 

with X's and F's real, and form linear combinations of the ordinary electromagnetic gauge 
potental with the eighth component of the colour gauge potential [14], 



V = A 8 cos6 + Asm9 

V = -A 8 sin0 + .Acos0, 



(2.20) 



where the "Weinberg angle" is given by 



2e 

tan# = — , (2.21) 

V3g 



t It is very important to maintain this general form of the CFL condensate, since its 
nontrivial winding onto the gauge group, SU(3) C x U(l) em gives rise to vortex filaments. 
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The Ginzburg-Landau theory and the surface energy 
Next we expand (2.11) to quadratic order in X's F's A : s and A. We find 
T = T min + J d s f{^(V x V) 2 + tr[(V x Wf + m 2 w W 2 } + ^[(Vx Z f + m 2 z Z 2 



+ l -[{VX) 2 +m 2 H X 2 ] + \[(VX l )(VX l ) + m' 2 X l X l ] + \{VY) 2 } 



(2.22) 



where W = Ya=i and 



W l = A l VY l for I = 1,2, ■■■7 

raw 

z = v - — vr 8 . 

m z 



(2.23) 



The masses of the excitations provide us with the relevant length scales which are the 
coherence lengths, £ and £' defined by: 

m 2 , = (6i + 36 2 )</> 2 = ^, m£ = &k/> 2 = p (2.24) 

that indicate the distances over which the di-quark condensate varies and the magnetic 
penetration depths, 5 and 5' by 

m| = 4sr 2 ( /> 2 1 sec 2 6> = , m 2 ^ = m|cos 2 = . (2.25) 

The excitations that violate (2.5) correspond to the Goldstone bosons associated with 
chiral symmetry breaking, and the r\ particle that becomes massive through the anomaly. 

In the rest of this section, we shall determine the critical magnetic field of a homo- 
geneous colour superconductor. For simplicity, the diamagnetic response of the quarks in 
the normal phase will be neglected. In the presence of an external magnetic field, the 
thermodynamic function to be minimized is the Gibbs free energy T, which is a Legendre 
transformation of the Helmholtz free energy T, i.e. 



r = r 



-H- J d 3 rV x A. (2.26) 

Following the decomposition (2.20), we write 

f = f 1 + f 2 (2.27) 

with 

d3 ^A/Z F ij F *J + 0^ X ^ 2 + 4t *(D®)hD$) +4atr$t$ + 5 ltr ($t$)2 

+ & 2 (tr$t$)2_£. Vx y sin #] (2-28) 
f 2 = J d 3 f[^(V xV) 2 - H -VxVcos6] 
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The minimization of r 2 gives rise to 

VxV = #cos#, (2.29) 

and no distinction will be made between the super-phase and the normal phase. In the 
remaining of the section we shall concentrate solely on r\. 

A homogeneous super-phase is characterized by a perfect Meissner effect, V x V = 
and the CFL condensate $ = (po, which gives rise to the minimum of T, F s = r m i n . given 
by (2.16). Fi for a homogeneous normal phase is minimized by V x V = Hsin9 and the 
minimum reads T n = — fi|if 2 sin 2 9. The critical magnetic field H = H c is determined 
from the equilibrium condition 

12o 2 

r n = r s = r min = -o (2.30), 

which leads to 



6(2 

Substituting the parameters (2.4) which were calculated using weak coupling approxima- 
tion into (2.30) yields 



"^^Tm^^-h^ (2 ' 32) 

Extrapolating this expression to a realistic value for the quark chemical potential, for 
example fx = 400MeV, and using the one-loop running coupling constant 

2 _ 127r2 

9 ~ CiN c -N f )In± (2 - 33) 
with N c = Nf = 3 and A = 200MeV, we find 

9 ^ -5.6° (2.34) 

and 

H c = 1.47 x 10 20 (Mf)(i_ Gauss, (2.35) 

which is much higher than the typical magnetic field inside a neutron star. This estimation 
of course depends on the validity of the extrapolation of the weak coupling formulas to 
realistic values of the chemical potential. 
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The Ginzburg- Landau theory and the surface energy 
3. The Domain Wall Problem. 



Let's consider an interface between the normal phase and the colour superconducting 
phase, for example the y — z plane. The equilibrium in the bulk is maintained by a uniform 
external magnetic field of critical magnitude, H — H c (. All quantities then in both phases 
will depend solely on x. It is evident from symmetry that the gauge field lies in one plane, 
let this be the x — y plane, so that V = V{x)y and V x V = The boundary conditions 
on the Ginzburg-Landau equations in the problem that we are considering (corresponding 
to the normal and colour superconducting phases as x i— > — oo and x i— > oo) are 



dV 

cp \— > 0, x l— 0) ~i — l—> H c s\n9 at x i— > — oo 
dx 

4> i— > 00) X l—> V^^o? -7- l— > at x I— > 00 



(3.1) 



The surface energy a is defined as 



a = Tl - Ts (3 2) 

Area of yz plane 

where the Gibbs free energy of a homogeneous superphase, T s is given by (2.30), and it is 
equal to that of the normal phase in the presence of H = H c . 

The minimization of f 1 generates a set of coupled nonlinear differential equations- 
the Ginzburg-Landau equations-subject to the boundary conditions (3.1). Now we divide 
the field variables in equation (2.8) into two groups, the first group contains $o 5 < ^ 8 7 ^. 8 7 A, 
while the second one includes A 1 with / = 1, 2, • • • 7. A closer inspection of the structure 
of reveals the absence of terms in the free energy (2.2) that are linear in the variables 
of the second group. Therefore a solution of the equations of motion exists in which only 
the fields of the first group acquire non-zero values. Since this ansatz implements the 
maximum symmetry allowed by the boundary conditions, we expect that it includes the 
solution that minimizes the free energy of the domain wall. 

By setting = A 1 = for / = 1, 2, ■ ■ • 7 and transforming the remaining variables 

we arrive at the following expression for the Ginzburg-Landau free energy functional (2.2) 

fx = J d 3 f ^(V x V) 2 + 4|(V - z-^sec^)^! 2 + 4|(V + i^=sec9V)x\ 

1 

2 



|2 



+ 4a(|0| 2 + | X | 2 ) + &i(H 4 + ^|x| 4 ) + b 2 (\<P\ 2 + | X | 2 ) 2 ) - H c ( ■ (V x V)sm6 



(3.4) 
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The condition for the colour-flavour locking ( that the 3x3 matrix $ is proportional to 
a unitary matrix ) amounts to \x\ 2 = 2|</>| 2 . Consequently, we shall use instead of the 
variable x and the functions V, (j> an d X the dimensionless quantities 

\J — 3ci 

s = Sx, <j> = 4>ou, \ = \pl(\)QV, V = Acos9. (3.5) 

The boundary conditions in terms of the new dimensionless quantities become 

u \— > 0, v i — ^ 0, A' \— > 1 at s i— > — oo 

(3.6) 

tt i— > 1, u i— > 1, 4 h0 at si-^oo 

where prime indicates differentiation with respect to s. The surface energy per unit area 
in terms of the dimensionless quantities becomes 



6a 2 
cr = —— 



J°°Js{\^ - I) 2 + ^(u 2 + 2</ 2 ) + \A\2u* + , 2 ) - \(u* + 2, 2 ) 



+ -3-(2n 4 + 2u 2 v 2 + 5v 4 ) + ±-p{u 2 - v 2 ) 2 } 

to to 



(3.7) 



where 

k= 5 - (3.8) 
is the Ginzburg-Landau parameter, b = \{b\ + 362), and 

h - 36 2 . , 

" = hTW 2 (39) 

is another dimensionless parameter. The weak coupling calculation fixes p to be equal to 
— 2- The Ginzburg-Landau equations which determine the profile of the colour condensate 
-u, v and the magnetic field in the colour superconductor are determined by minimizing 
the surface energy with respect to functions w, v and A 

- A" + -A(2u 2 + v 2 ) = 0, 



— u" + (A 2 - l)u + -(2u 2 + v 2 )u + -p(u 2 - v 2 )u = 0, (3.10) 
k 3 3 

1// l,,n , s I/O _On 1 



^rv" + -(A 2 - A)v + -(u 2 + 5v 2 )v - -p(u 2 - v 2 )v = 0. 
k 4 6 6 

In the presence of a non-zero A, as it is required by the boundary conditions, this set of 
equations does not admit a solution in which u = v everywhere. Stated differently, in the 
presence of an inhomogeneity and a non-zero gauge potential the unlocked condensate- 
the octet-has to acquire a non-zero value somewhere. This situation is analogous to the 
Ginzburg-Landau theory of a cuprate superconductor [15], where the s-wave condensate 
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becomes non-zero in the vicinity of a vortex filament while the <i-wave condensate domi- 
nates in the bulk. It is easily verified that equations (3.10) have the first integral 

l -A' 2 + (u 2 + 2v' 2 ) - l -A\2u 2 + v 2 ) + \(u 2 + 2v 2 ) - ^(2n 4 + 2u 2 v 2 + 5v*) 

1 1 ' 

- T8 p(u -v) =- 

(3.11) 

which implies, according to the boundary conditions (3.6), that 

A(oo) = 0. (3.12) 



Let cr m in be the minimum surface energy. It is obtained from the substitution of the 
field variables u, v, A that satisfy the GL equations (3.10) into (3.7), and is a function of 
the parameters k and p. We proceed now to prove the following lemma 

Lemma: 

(^). £ * (3 - 14) 



Proof: 

According to (3.7), 
da ■ \ 4a 2 r °° 



Ok ) p bn 3 



p 

oo 

+ 



ds(u' 2 + 2v' 2 ) 

f°° ^ r <fomin t dA(s) \ So-min f dll(s) \ &X mi n / 8v{s) \ j ^ ^ 

J_ 00 IdA(s) V 8k J p 5u(s) \ 8k J p 8v(s) V 8k ) p . 

The second term vanishes because of the equation of motion (3.10) and the boundary 
conditions (3.6), (3.12). Consequently, we find that 

(^X=-S/_>"' 2+2 ^°- < 3 - i6 > 

Similarly 

(^). = < 317 » 

The lemma is then proved. 
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Furthermore lets denote the critical value of the Ginzburg-Landau parameter for which 
Omin vanishes by K c (p). We shall prove then the following theorem 

Theorem: 

«c(p) < ^= (3-18) 

and 

^ > 0. (3.19) 
dp 



Proof: 

Consider a special field configuration u = v, which satisfies the boundary conditions 
(3.6) but not the Ginzburg-Landau equations (3.10) with 0. Therefore 

<W < ^ f°° ds \\{A> - l) 2 + \u' 2 + l -A 2 u 2 -u 2 + \u 4 } = a (3.20) 

J-oo Lz K I Z J 

The minimization of a with respect to u and A yields the set of differential equations, 

- A" + Vu 2 = 0, 

1 1 (3-21) 

-\u" + Ua 2 -2)u + u* = ^ 

K Z 

subject to the boundary conditions (3.6). We recognize that (3.20) and (3.21) correspond 
exactly to the domain wall problem for an ordinary superconductor analyzed by Ginzburg 
and Landau in their original work. It follows then that a = at k = A=. On the other 
hand, this field configuration does not satisfy the equations of motion (3.7) and therefore 
Omin < at k = Following (3.13) of the lemma we establish the first equation of the 
theorem. 



The second equation of the theorem follows from the identity 




(3.22) 



and the lemma. The theorem is then proved. 

The solution to the Ginzburg-Landau equations (3.10) and the critical value of the 
Ginzburg-Landau parameter k for various values of p G (— l,oo) were also determined 
numerically. The continuous range s G (—00,00) was replaced by a finite lattice with 

s n = —L + ne (3.23) 
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where n = 0, 1, 2, N+l and (N + l)e = 2L. A good approximation amounts to e and L 
that satisfy e << min(l, 1/k) and L » max(l, The field variables u n , v n , and A n 

are assigned to each site and their derivatives can be approximated by u' n = (u n+ i —u n )/e, 
v 'n = ( v n+i — v n )/e, and A' n = (A n+ i — A n )/e. The surface energy is then written 



a 2 N i i 

- = ^ E & A 'n ~ ^ + ^2« 2 + 2 < 2 ). 



n=0 

AT+1 111 1 

+ E e - [g^ 2 ^ + ^ - 3 ( w n + 2t£) + ^(2< + 2«Jt;; + 5t£) + _ ^2 )2 j } 

n=0 

(3.24) 

with eo = e7v_|_i = |e and e n = e for n = 1, • • • N. The multivariable function (3.24) is then 
minimized by iteration subject to the conditions u$ = v o = 0, un + i = vn + i = 1, A' Q = 1 
and Atv+i = 0. 

The critical GL parameter, k c versus p is plotted in Fig.l. The critical value of 
GL parameter in the weak coupling approximation (p = — |), is k c = 0.589. The two 
components of the di-quark condensate, the colour-flavour locked and the colour-flavour 
octet, and the corresponding magnetic field, that correspond to the solution with p = — | 
and plotted in Fig. 2. The numerical results are consistent with the inequalities 

(3.18) and (3.19) that we derived analytically. 

Substituting the weak coupling expression of the parameters, the coherence length 
reads 

2 _7C(3)/ 1 Y T c 

and the penetration depth 

: ' 7r /; cos 2 # (3.26) 



2a s p 2 T c -T 

2 

with a s = The Ginzburg-Landau parameter is given then by 



S 2 / 72tt 3 k B T c n 16.3 k B T c 
K= a = \hrl^ — cos# = ~=-^- 3.27 

If we consider a realistic value for the chemical potential, for example p = 400MeV, it 
follows from (1.3) that the colour superconductor will be of type I if ksT c < 14MeV and 
of type II otherwise. Furthermore if we extrapolate the weak coupling formula for T c [16], 
[17], [18], [19], [20], which is valid at asymptotic densities, to p = 400MeV we find ksT c = 
3.5MeV and the colour superconductor is of type I. Of course if nonperturbative effects 
raise T c significantly, for example by one order of magnitude, the colour superconductor 
could be of type II. 
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4. Concluding Remarks. 



In this paper we have applied the Ginzburg-Landau theory to the problem of a do- 
main wall which separates the normal phase from the colour superconducting phase of 
dense quark matter. The main purpose of this work was to establish a criterion that 
determines whether the colour superconductor is type I or type II with respect to an ex- 
ternal magnetic field. We initially derived the Ginzburg-Landau equations of motion for 
the gauge fields and the order parameter by minimizing the surface energy of the interface 
between the two phases. The condition that the surface energy is a minimum is equivalent 
to the requirement that the normal and the colour superconducting phases are in stable 
equilibrium with each other. This equilibrium is maintained by an external magnetic field. 
Using symmetry arguments we were able to simplify the original set of equations and 
rewrite them in terms of only two components of the order parameter and one gauge po- 
tential. The equations also contain two dimensionless parameters: the Ginzburg-Landau 
parameter k and p. Therefore the problem is more complicated than the problem of the 
metallic superconductor which was addressed by Ginzburg and Landau. Although we were 
unable to derive an analytical result analogous to the Ginzburg-Landau criterion (1.1), we 
derived several rigorous inequalities about the critical parameter k c , the value for which 
the surface energy vanishes, and its derivative with respect to p. Those analytical results 
are supported by a numerical solution of the Ginzburg-Landau equations. By extrapolat- 
ing the weak coupling approximation formulas for the parameters to realistic values of the 
chemical potential, for example fi = 400MeV, we find that the colour superconductor is of 
type I (type II) if the transition temperature is below(above) 14MeV. 

The value of our work is mainly theoretical. It addresses the broken U(l) gauge 
field onto which the electromagnetic field has a small projection. Although we have been 
refering to the typical chemical potential of a neutron star, the implications of our work 
might be quite limited [14]. First of all, there has not been any evidence yet that colour 
superconductivity is realised in the core of a neutron star. Even if this is the case, the 
critical magnetic field (2.35) is too high for the distintion between type I and type II to 
be observed unless the pairing force is strong enough to reduce the lower critical field by 
several orders of magnitude. In addition the small mixing angle (2.34) and the dependence 
of the magnetic response on the thickness of the crust of the CSC core make it difficult to 
observe the partial magnetic Meissner effect and the appearance of vortex filaments. 

Because of the strength of the QCD coupling and the ultra-relativistic Fermi sea, the 
fluctuations of the gauge field could be more significant than those in a non-relativistic 
superconductor [21]. These fluctuations contribute to the Ginzburg-Landau free energy 
functional an energy density 

Ar ~ k B T c 5~ 3 . (4.1) 
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As the condensate energy (2.16) is proportional to (1 — jr) 2 while 5 



1 — Tjr-, the 



contribution from the fluctuations will eventually dominate as T approaches T c and modify 
the nature of the phase transition. 

The fluctuation energy is estimated in the appendix within the framework of the 
Ginzburg-Landau theory, where we demonstrate that the inclusion of this energy modifies 
the phase transition from second order to first order. The parameter 



which describes the deviation of the first order transition temperature T from the second 
order one T c , is used as a measure of the importance of the fluctuations. We find that for 
realistic values of the quark chemical potential, fx = 400MeV and /csT c > 13MeV, t < 0.1. 
In that case, the Ginzburg-Landau free energy functional (2.2) represents a reasonable 
approximation. 

For t ~ 1 the approximation employed to derive the Ginzburg-Landau free energy 
(2.2) with (2.4) from QCD or from the Nambu-Jona-Lasinio effective action breaks down. 
Higher powers of the order parameter have to be included. It would be very interesting to 
develop a systematic weak coupling approximation in this case and we hope to be able to 
report progress in the future. 
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6. Appendix. 



In this appendix we shall estimate the effect of the fluctuations of the gauge field. A 
systematic evaluation of the contribution from the fluctuations to the free energy requires 
a resummation of a set of ring diagrams of QCD at finite temperature. Here we shall 
follow the treatment of Bailin and Love [21] and estimate the fluctuation terms within the 
framework of the Ginzburg-Landau theory. Due to the Debye screening the contributions of 



t 




(4.2) 



15 



The Ginzburg- Landau theory and the surface energy 



the longitudinal components of the gauge potentials are supressed, and only the transeverse 
components contribute in the static limit. We expect that our simple-minded estimation 
captures the fluctuation terms to the leading order of the condensate. 



Assuming that the small fluctuations are governed by the quadratic form of the 
Ginzburg-Landau free energy functional, the shift of the free energy due to the Meiss- 
ner effect is given by 

Ar = -k B T In J rr L JTT J r L 7 J g 6.1 
f[dW][dZ]e- s ° ' 

— * — * 

where W and Z represent the transverse degrees of freedom and 



S = | J rf 3 f[(V x W) 2 + m 2 w W 2 + (V x Z) 2 + m\Z 2 \ 
So = ^J d 3 f\(V x W) 2 + (V x Z)\ 
Completing the gaussian integral, we find that 

- kBTc (7S'- 3 + 5~ 3 ) = -^(7 + sec 2 ^)^T C |^| 3 . 



(6.2) 



(6.3) 



67T 37T 



In the previous calculation, dimensional regularization was employed to eliminate ultravi- 
olet divergences. By including this term to the homogeneous Ginzburg-Landau free energy 
with $ = </>, we find 

4o 3 

T = 12a\(j)\ 2 + 12b\(j)\ 4 - — (7 + sec 2 6)k B T c \(j)\ 3 . (6.4) 

37T 

The presence of the cubic term will induce a first-order phase transition at T > T c The 
temperature is determined by the condition that the free energy curve T vs. \4>\ becomes 
tangent to the <^>-axis at some ^ 0, i.e. 

a= 1 -^(7 + sec 2 e) 2 (k B T c ) 2 . (6.5) 

Substituting the weak coupling expression for a and b we obtain that 

t _ T 49C 2 (3)(7 + sec 2 ^) 2 « 3 » 2 

T c 139968tt 5 k B T c [ ' ' 

With \i = 400MeV, and k B T c > 13MeV, t < 0.1. As a result the fluctuations are 
not signigicant and the Ginzburg-Landau energy functional (2.2) represents a reasonable 
approximation. 
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Figure 1. The critical Ginzburg-Landau parameter k c vs. the parameter p. 




Figure 2. The solutions to the Ginzburg-Landau equations (3.10) for p = — | and k = 0.589. 
The solid line represents the colour-flavour locked component, \{u -\- 2v), the dashed line 
represents the colour-flavour octet component, — it) and the dotted line the colour- 

magnetic field, A'. 
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